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Abstract 



We show that the angular distributions of leptons or jets due to it spin 

correlations allow a determination of the top chromomagnetic moment k with 

an accuracy of order 0.1. The method is quite insensitive to background, event 

reconstruction, and other experimental uncertainties. The total event number 

\Q ' is important, so we suggest the inclusion of all-hadronic events. 

o" 
°\ 
Sj ; 1 Introduction 

A primary goal of high energy experiments is to observe deviations from the standard 

model. Prior to the detection of, or in the absence of, new light resonances such as 

a Higgs particle, one is left with low energy remnants of loop effects, such as higher 
• i-^ ■ 

/\ dimensional operators. The dimension-5 top chromomagnetic moment operator is of 

H ' 

special interest due to its relationship with the top massQ 

In dynamical schemes of the top mass the estimate is k ~ m 2 /A 2 , where A character- 
izes some scale of physics from which the top is deriving its mass. This is typically 
the electroweak symmetry breaking scale or smaller. The assumption here is that 
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the heavier degrees of freedom also carry color, in which case the chromomagnetic 
moment is related to the top mass diagram by attaching a gluon. The situation is 
different for the standard model or other models involving elementary scalar fields. 
In these models the chromomagnetic moment arises at one loop, in contrast to the 
tree level origin of the top mass, and this leads to a 1/(Att) 2 loop factor suppression 
of the chromomagnetic moment. 

The effect of a top chromomagnetic moment in hadron collider experiments has 
been investigated |l|. Naively, one would expect that the energy dependence of 
the new contribution differs significantly from the standard model contribution, as 
the derivative acting on the top field introduces a momentum factor. Actually the 
structure of the contributions is schematically like 1 + k + k 2 - 2 -?. The would-be 
enhancement factor -2- in the k term is absent, since the interference between the 
standard chirality conserving amplitude and the chirality changing (k) amplitude 
requires an additional m t factor. This, together with the fact that most it pairs are 
produced not far above the threshold, explains the fact [|IJ that various differential 
cross sections in top momentum variables are sensitive to n essentially only through 
an overall change in their normalization. The extraction of k from such observables 
at the Tevatron is thus difficult because (i) the luminosity is not known with a good 
accuracy, and (ii) other operators may also contribute to the structureless excess. 

We are led consider other observables. In particular a chirality changing contri- 
bution should affect the spin correlations of the it state, thus modifying the angular 
distributions of the top decay products. Analyses of the it spin structure in the dom- 
inant qq — > g — > it process in the standard model are provided in Refs. []2], 0, f|, [|. 
Some of the top spin information resides in the direction of the charged lepton (or 
down-type quark), the decay product of the top in the chain t — > W — > I. When only 
this direction is measured in the top rest frame, relative to the top spin direction s, 
the spin correlation shows up as [|, [| 
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d is the angle between s and the charged lepton direction, and similarly for $ and 
s on the t side. (8 is the top scattering angle.) Nn (N x ) is the number like-spin 
(unlike-spin) it pairs (by like-spins we mean spins both parallel or both antiparallel 



to their respective axes). N\\ vanishes for a particular choice of s, which defines the 
spin basis dubbed "off-diagonal" 0. Then the t and 1 spins are 100% correlated, up 
to small corrections from the gluon induced process. 

Unfortunately the $ and d variables as defined in the off-diagonal basis provide a 
poor probe of the chromomagnetic moment operator. We have calculated the differ- 
ential cross section for qq —>■ g — > it — > b UD + bUD to O(k) in terms of an enlarged 
set of phase space variables. (At Tevatron energies 0(n 2 ) effects may safely be ig- 
nored.) Requiring that both the top quarks and the W ± be on shell we find, again 
with variables defined in the off-diagonal basis, 

<I<J <xK K =KW+KMK (3) 
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(p and ^ are the azimuthal angles around the corresponding spin axes s and s and 
7 = (1 — (3 2 ) where (3 is the top velocity. 

We immediately observe that an integration over the azimuthal variables gives 
back 
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This result clearly elucidates the relationship of our approach to both Ref. |1| and ||. 

The first factor contains all information used in [Q], with the 1/2 factor suppressing 

the dependence on both f3 and 9. The second factor is the same as found in [0 

(Eq. (H) above in the off-diagonal basis), so that the i? and ■# variables by themselves 

introduce no dependence on k. 

The complicated structure of the cross section formula prompts us to look for a 
simple variable which still carries most of the information useful for the extraction 
of k from a data sample. Attempts to find such variables were discussed in M. But 
the structure of Eq. (H) is such that all useful information is encoded in the sin's and 
cos's of the azimuthal angles which change sign as we move around the phase space. 
The best variable should incorporate as much of this information as possible. 

Although not necessarily simple, this "best" variable is provided by the maximum 
likelihood method, which is the standard way of extracting the value of a parameter 



of a distribution from a sample. Noting that 1Z K (Eq. (|3|)) is defined to be linear in 



n, this variable is 



The maximum likelihood method uses only the values of this variable in the determi- 
nation of K. 

In order to assess the observability of k in Run 2 of the Tevatron we built a 
Monte-Carlo generator and used a maximum likelihood estimate on the generated 
data sample to determine the possible statistical accuracy of the estimated value of 
K. We found, taking a reasonable 1500 events, a "1-sigma" accuracy of O(±0.1). 
What is quite remarkable is that the accuracy of the determination is very stable 
against factors like (i) the introduction of various cuts, (ii) crude inaccuracies in the 
measurement of the parameters of each individual jet, (iii) the inability to tell apart 
the two hadronic decay products of the W or (iv) the introduction of a background 
even as large as ten times the it signal. All these insensitivities are due to the 
complicated but smoothly varying function in Eq. (||), which is not easily faked or 
influenced by any of the above factors. 

The only factor of any importance is the number of it events. We find that a 
sensible determination of k requires 1000-2000 events and any increase in the event 
number significantly improves the determination. However, for each individual event 
all that is needed is a very crude estimate of the directions of all six decay products 
(on the partonic level), which is even possible in the case of all-hadronic decay. It 
seems advantageous to relax the standard set of cuts, at the cost of larger backgrounds 
and loss of precision in each event, in order to gain in the event number. The method 
may then provide an estimate of k superior to the determination from the total cross 
section. In any case, the method advocated here gives a completely independent way 
to help determine the source of any observed cross section anomaly. 

The rest of the paper is organized as follows. Sec. |2| describes our calculation of 
the spin density matrix, explains how we handle the six-particle phase space, and 
produces our result for the differential cross section. Sec. ||| describes the maximum 
likelihood estimate method while gives Sec. ^ gives the numerical results. 



2 The structure of the differential cross section 



At the Tevatron the dominating process for tt production is quark fusion qq 
tt. The amplitude of this process at tree level is 
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which directly translates into a spin density matrix for the top 
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where color summation/averaging is understood, as well as averaging for incoming 
quark spins. The indices s, s refer to the top (antitop) spin, and s', s' are the corre- 
sponding spin indices in the conjugate amplitude. They are all measured with respect 
to an arbitrary (for the moment) spin direction s M , s M . For defmiteness, we choose the 
spatial part of the s M , s M vectors back-to-back at an angle ip to the beam direction in 
the zero momentum frame (ZMF) frame (see Fig. [1]). 




Figure 1: The definition of a spin basis in the ZMF frame. The z-axis points out of 
the page and is used to define the azimuthal angles <p and Tp. 

We calculate each of the elements of the density matrix. Its structure is greatly 
simplified in the "off-diagonal" basis defined by || 

tan-0 = — — r with < ib < — . (8) 



In this basis the standard model contributions from like-spin tt pairs vanish. We 
determine how the result in || is modified at linear order in k: 
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We use the shorthands sg = sin.9 and eg = cos 9. The T states are those where the 
t (t) have spin projection +1/2 on their respective spin axes. Note that this matrix 
is rank 2, as in the standard model case, so it involves only two tt spin states. Any 
other choice of the angle ip (such as if) = for the helicity basis) only makes the form 
of p more complicated. 

All the K-dependence of the process is encoded in this density matrix. Note that 
Trp = 4(1 + 2k — ^(3 2 Sg) is nothing but the first factor in Eq. (f|). The information 
in the off-diagonal elements will be observed in the angular distributions of the decay 
products. 

The correct choice of the variables in the six particle phase space (see Fig. 0) should 
keep the expressions for both the amplitude and the six-particle phase space measure 
as simple as possible. We should also make sure that the regions of integration are 
simple. For example, no combination of angles and/or energies in the ZMF frame 
will satisfy these requirements. For these reasons we have chosen the following set of 
variables, in addition to the ones that define the ZMF frame: 

• The scattering angle 9 between q and t in ZMF and the corresponding (irrele- 
vant) azimuthal angle defining the scattering plane. 

• The angle d in the top frame between the spin direction s and the D quark 
(or charged lepton) and the corresponding azimuthal angle (p. Note that tp — 
corresponds to the positive z direction in Fig. [1|. 

• The angle $ in the antitop frame between the spin direction s and the D quark 
(or charged lepton) and the corresponding azimuthal angle Tp. Note that Tp = 
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U(v) 




D(/ ) 



U(v) 



Figure 2: The full cascade of tt production and decay. The decay products of the W + 
are an up-type U = u,c quark (or a neutrino) and a down-type antiquark D = d,s 
(or a positively charged lepton, e + or /i + ). 

also corresponds to the same positive z direction in Fig. |l|. 

• The angle A in the top frame between the direction of W + and the D quark 
and the corresponding azimuthal angle 0; also the respective angles A, <fi on the 
t side. 

These ten variables completely describe the six-particle phase space (which is 
10 dimensional due to four restrictions from momentum conservation and four from 
the on-shell conditions of the t,t, W^). Their corresponding ranges are completely 
independent of each other and the phase space measure is very simple. We emphasize 
though that these angles are not all defined in the same frame. 
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We calculated the decay rate of each of the 16 spin states of the tt system indi- 
vidually and used a Breit-Wigner formula to find the qq — > 6 quarks cross section. 
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Here /?w = (m| — m^)/(mf + 77t|f) and we note how the integration over the four 
variables A, 0, A, can be trivially done. We find 
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where the two matrices represent the angular dependences of the t and t decay rates. 
The structure of the above matrices implies that the off-diagonal elements of p gets 
encoded in the p and ^-distributions. By integrating out p and Tp we would see only 
the diagonal elements of p. Substituting Eq. @ into Eq. ( pTTj) gives the result for TZ K 
in Eq. (|3]). Thus we have managed to express all spin correlation information in terms 
of five angles and the top quark velocity. 

3 Maximum likelihood and the best variable 

We would like now to assess the statistical errors in the estimation of k at Run 2 of 
the Tevatron. For that purpose we simulate the actual process by a sample of events 
generated by a "theorist's Monte- Carlo" , an event generator which uses the quark- 
level differential cross sections in Eqs. (^,|10|) folded with MRST parton distribution 
functions. For our purposes of crude estimation it is sufficient to treat gluon induced 
processes as background and to ignore higher loop and radiative corrections. 

As the next step we took the generated sample with a given small k = n true and 
used an unbinned maximum likelihood method to estimate k from the sample. We 
calculate the log-likelihood function 

log£(«:) = ]riogp« (12) 



where p^' is the probability for an event to occur at the point of the phase space 
labeled by i. This probability obviously depends on the measure used in the phase 
space but one observes that all the ambiguity in redefining the phase space variables 
(and so, the measure) and in redefining the grid used to produce a nonzero probability 
for each event produces only a K-independent Jacobian factor multiplying C(k). We 
remove this ambiguity altogether by subtracting the standard model value \ogC(n = 
0) from Eq. (|12|)- This will not affect the shape of the likelihood function nor the 
position of its maximum. 

The quantity p K , however, involves a nontrivial normalization issue. As a proba- 
bility it should be normalized to one: 

(l) f_m+Kfm _ m) l +^iM (u) 

Recall that we are considering the amplitude-squared to linear order in k. It is a 
function of the phase space variables collectively denoted as $. The first factor in 
Eq. (|13D, being independent of k, can be absorbed into the normalization ambiguity 
of the likelihood function, while the last denominator cannot. It is the reflection of 
the fact that we estimate k on a sample of given size, so we loose the information 
related to the total rate. On the other hand we are keeping the information residing 
in the angular distributions, and in this way our approach is orthogonal to that in M. 
The probability of each event in Eq. ( O ) is expressed in terms of one variable 
F (®) = j§y and a constant A = j^J^ , so 

log£( K )=Elog^f. (14) 

The variable F carries all information used by the maximum likelihood method and 
is thus the "best" variable appearing in Eq. (|5|). 

4 The stability of the estimation 

We estimated k in samples of various sizes, with various true values of k fed into the 
Monte-Carlo generator. As can be understood from Eq. ([14]) , the likelihood function 
becomes singular at K — —F^ 1 for each event i. These singularities are an artifact 
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of our linear approximation in k, and tell us that our estimation is valid only when 
K <C min |-Fj _1 |. But this turns out to cover the theoretically interesting range of k. 




Figure 3: The likelihood function for a specific Monte-Carlo run with n true = 0, based 
on 1500 events. The horizontal line represents the 1-sigma interval. The skewing of 
the function is common to the different runs, but the location of the maximum is of 
course randomly distributed about k — 0. 

Fig. |3| shows the likelihood function for a characteristic example run with 1500 
observed it events. Following tradition, we will quote the interval where the log- 
likelihood function drops by 1/2 as a 1-sigma interval, although we have to keep in 
mind that such a statement is mathematically inaccurate. This interval is shown as 
the horizontal line. We see that the estimation in Fig. |3| is consistent with the true 
value of k and corresponds to a 85% probability level. 

We ran the Monte-Carlo with various different initial settings of the random gen- 
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erator. We found the following 1-sigma intervals by averaging the values: 

iV = 200 Ak = 
iV = 900 Ak = 
iV = 1500 Ak = 
N = 2000 Ak = 
N = 5000 Ak = 



+0.30 
-0.15 

+0.22 
-0.12 

+0.10 
-0.07 

+0.10 
-0.07 



(15) 



+0.05 
-0.03 



We immediately observe the importance of the event number. The number of 
observed ti events on Run 2 is expected to be (9(1500). With that number the 
accuracy is in the order of An = 0.1, not very far from the prediction of some models 
of new physics. 

We next address the question whether experimental realities significantly deteri- 
orate these results. We considered the following issues: 

• In the hadronic decay mode of the W it is almost impossible to tell which of 
the two jets comes from which quark. We set up a similar maximum likelihood 
estimation of k where the new probability function now accounts for a 50% 
probability of misidentifying the two jets. The change from the values in ([15]) 
is insignificant (by this we will mean a less than 20% change in Ak.) This is 
true for "symmetrization" on only one side (i.e. for the W + only, relevant for 
semileptonic events) and for both sides, relevant for all-hadronic decays. 



We imposed several sets of cuts on the data to see their effect on the extraction. 
The variable Fi is unmodified but the normalization constant A = J r ,,,,., 
is modified as the integration region is changed. We found a value for A by 
estimating k on a very large sample and requiring that it reproduce K = K true . 
For example, we imposed the same parton level cuts that had been used in the 
analysis of the semileptonic events in Run 1 |§ , (a charged lepton with E? > 
8GeV, \rj\ < 1.0, a neutrino with E T > 20GeV, three jets with E T > !5GeV 
and \rj\ < 2.0 and one more jet with Et > 8 GeV and \r]\ < 2.4. We found only 
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insignificant changes in the accuracy when the number of events passing the 
cuts were fixed. Imposing a jet separation cut of Aa > 20° similarly had little 
effect. 

There are usually very large uncertainties in extracting the quark three mo- 
menta. We set up a random distortion of the energy and direction of each jet 
by a relative portion of 5 = 30%. Imposing the top and W mass conditions on 
the data we found again only insignificant changes. As a check, S = 90% did 
destroy the usefulness of the data as would be expected. 

Finally, we included a very crude model of the background, for which we used 
a uniform distribution in the angular variables. For the number of it events 
constant, iV = 1500, the size of the 1-sigma interval increases as follows: 



N< 



total 



N< 



total 



15000 An 



4500 An 



Ntotai = 3000 Ak 



N i 



total 



Ni 



total 



1875 



1500 



An 
An 



+0.15 
-0.09 

+0.13 
-0.08 

+0.11 
-0.07 

+0.11 
-0.07 

+0.10 
-0.07 



(16) 



This slight increase, however, shows only that even a large background can be 
tolerated in exchange for a larger event number. 

We see that many of the experimental difficulties have little effect in our case, 
basically due to the fact they cannot easily mimic the intricacies of the angular 
distribution in Eq. @. Systematic errors would similarly be dangerous only if they 
include a component with a correlation between the various angles similar to Eq. fl5|). 

In conclusion, we have shown how it spin correlation information can be extracted 
from the angular distributions of the it decay products. This provides an independent 
and competitive way to observe new physics contributions to the chromomagnetic 
moment of the top quark. The statistical limitation on the extraction of k at Run 2 of 
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the Tevatron is (9(±0.1), mainly coming from statistics. The extraction is surprisingly 
stable against experimental uncertainties and we suggest that cuts should be relaxed 
as much as possible in order to gain in the event number. The inclusion of all-hadronic 
events should also be given serious consideration. 
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